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We advocate a set of approximations for studying the finite temperature be- 
^. havior of strongly-coupled theories in 0+1 dimensions. The approximation 

r-| ! consists of expanding about a Gaussian action, with the width of the Gaus- 

sian determined by a set of gap equations. The approximation can be applied 
]S^ [ to supersymmetric systems, provided that the gap equations are formulated 

j_] ■ in superspace. It can be applied to large-iV theories, by keeping just the 

planar contribution to the gap equations. 

We analyze several models of scalar supersymmetric quantum mechanics, 
and show that the Gaussian approximation correctly distinguishes between 
a moduli space, mass gap, and supersjTiimetry breaking at strong coupling. 
Then we apply the approximation to a bosonic large-A^ gauge theory, and 
argue that a Gross- Witten transition separates the weak-coupling and strong- 
coupling regimes. A similar transition should occur in a generic large-A^ 
gauge theory, in particular in 0-brane quantum mechanics. 



1 Introduction 



Recent developments have made it clear that at a non-perturbative level 
string or M-theory often has a dual description in terms of large-N gauge 
theory [|I], ^. This has many interesting implications, both for string theory 
and for the behavior of large-N gauge theory. These developments have been 
reviewed in 0. 

A general property of these dualities is that semiclassical gravity can 
be used only in regimes where the gauge theory is strongly coupled. For 
certain terms in the effective action, which are protected by supersymmetry 
0, a weak coupling calculation in the gauge theory can be extrapolated to 
strong coupling and compared to supergravity |^, ||, 0; H- But in general 
non-perturbative methods must be used to analyze the gauge theory in the 
supergravity regime. Although this reflects the power of duality, as providing 
a solution to strongly coupled gauge theory, it is also a source of frustration. 
For example, one would like to use the duality to understand black hole 
dynamics in terms of gauge theory @, |10|, 0, |l|, |l|, |l|, |15|, [l|, [l^, [T|, 



lOj pot pi| , P^ . But progress along these lines has been hampered by our 



inadequate understanding of gauge theory. 

Thus we would like to have direct control over the gauge theory at strong 
coupling. Such control is presumably easiest to achieve in a low-dimensional 
setting. The specific system we have in mind is 0-brane quantum mechanics 



23 1 , in the temperature regime where it is dual to a 10-dimensional non- 



extremal black hole |p4[. The dual supergravity makes definite predictions 
for the partition function of this gauge theory. In the temperature regime 
under consideration, the partition function is in principle given by the sum of 
planar diagrams. Although directly summing all planar diagrams seems like a 
hopeless task, one might hope for a reasonable approximation scheme, which 
can be used to resum a sufficiently large class of diagrams to see agreement 
with supergravity. 

In this paper we develop a set of approximations which, we hope, can 
eventually be used to do this. Essentially we self-consistently resum an infi- 
nite subset of the perturbation series. These methods have a long history in 
many-body physics ||2^ and field theory [^, ^, and have been used to study 



QCD [^, ^, large-N gauge theory |^, and even supersymmetry breaking 



5T| . The basic idea is to approximate the theory of interest with a Gaussian 
action. The Gaussian action is (in general) non-local in time, with a variance 
that is determined by solving a set of gap equations.^ One can systematically 
compute corrections to this approximation, in an expansion about the Gaus- 

^We refer to gap equations even when studying models without a mass gap. 



sian action. To study large-N theories in this formahsm, one keeps only the 
planar contribution to the gap equations. To study supersymmetric systems, 
we will see that one needs to formulate the gap equations in superspace: loop 
corrections to the auxiliary fields must be taken into account in order for the 
approximation to respect supersymmetry. 

An outline of this paper is as follows. In section 2 we illustrate the ap- 
proximation in the simple context of + dimensional (p^ theory, and discuss 
various ways of formulating gap equations. We also discuss (but do not re- 
solve) the difficulties with implementing gauge symmetry in the Gaussian 
approximation. In section 3 we turn to a series of quantum mechanics prob- 
lems with A/" = 2 supersymmetry but no gauge symmetry, and show that 
the Gaussian approximation captures the correct qualitative strong-couphng 
behavior present in these simple systems. In section 4 we discuss large-A" 
gauge theories in -|- 1 dimensions, and argue that generically a Gross- Witten 
transition is present, which separates the weak-coupling and strong-coupling 
phases. In section 5 we apply the Gaussian approximation to two large-A^ 
quantum mechanics problems: a bosonic gauge theory, and a supersymmetric 
matrix model with a moduli space. We discuss the relevance of these results 
for 0-brane quantum mechanics. 



2 Formulating the Gaussian approximation 



2.1 A simple example: 0^ theory in + dimensions 

Our approach to studying strongly-coupled low-dimensional systems can be 
illustrated in the following simple context. Consider + dimensional 0^ 
theory, with action 

1 /I ,. 1 



The partition function is given in terms of a Bessel function, by 

Z = e-^^ = r d<j>e-' = -^ey^'^Ky, f-1- ) . (2) 



V2 "^ \Sg 

Expanding this result for weak coupling one obtains an asymptotic series for 
the free energy. 

PF = -i log (2vr/) + |/ - 3/ + ^/ + 0{g') (3) 

This series can of course be reproduced by doing conventional perturbation 
theory in the coupling g'^. 



But suppose we are interested in the behavior of the free energy at strong 
couphng. From the exact expression we know that this is given by 

PF = -\ \ogg - log ^^^ + 0{l/g) . (4) 

Is there some approximation scheme that reproduces this resuh? Clearly 
working to any finite order in conventional perturbation theory is hopeless. 
But it turns out that by resumming a subset of the perturbation series one 
can do much better. To this end consider approximating the theory (P with 
the following Gaussian theory. 



S. = ^y (5) 

The width of the Gaussian o"^ is for the moment left arbitrary. By writing 

one obtains the identity 

13F = (3Fo -< e-^^-^"^ - 1 >c,o (6) 

where the subscript C, denotes a connected expectation value in the Gaus- 
sian theory (j^). By expanding this identity in powers oi S — So we obtain a 
reorganized perturbation serie^ 

PF = pFo+ <S-So>o-l<iS- SoY >c,o + ■ ■ • • (7) 

The first few terms are given by 
PFo = --log2na^ 
<S-So>=-{—-l] + 



2 V^^ / 4^^ 

The identity @ holds for any choice of 5*0 . But we wish to truncate the 
series (0), and the truncated series depends on the value of a^. So we need a 

^ There is a real advantage to keeping at least the first two terms in the reorganized 
perturbation series: one is free to add a constant term to Sq, which may depend on the 
couplings or temperature but is independent of the fields. Any such additive ambiguity 
cancels between BFq and < S ~ Sq >o. 



+ 



Q 



Figure 1: The gap equation for the cf) propagator. Heavy hues are the dressed 
propagator o"^ and thin hues are the bare propagator g"^. 

prescription for fixing cr^. One reasonable prescription is as follows. Consider 
the Schwinger-Dyson equation^ 

This is an exact relation among correlation functions in the full theory. It 
seems reasonable to ask that the same relation hold in the Gaussian theory, 
that is, to require that < ^ {cj? + 0^) >o= 1. This implies that cr^ obeys the 
gap equation 

a^ g^ g^ 

This equation has a simple diagrammatic interpretation, shown in Fig. 1: it 
self-consistently resums certain self-energy corrections to the propagator. 

At weak coupling the gap equation implies o"^ = (7^ + 0{g^) and thus the 
leading term in the Gaussian approximation is 

PFo = -hog2ng' + 0{g^). 

This matches the weak coupling behavior seen in (|^). At weak coupling 
the higher order terms in the reorganized perturbation series are small, sup- 
pressed by increasing powers of the coupling. 

Of course it is not surprising that the Gaussian approximation works 
well at weak coupling, since at weak coupling the action ([^) is approximately 
Gaussian. But what is rather remarkable is that the Gaussian approximation 
also gives good results at strong coupling, where the 0^ term in S dominates. 
At strong coupling the gap equation gives ^^ = ^ + ^(1) ^^^ thus the 
leading term in the Gaussian approximation is 

PFo = -Uogg ~^\og^ + Oil/g) . 

Thus i3Fq has the same leading strong-coupling behavior as the exact result 
(H). Moreover, by reorganizing the perturbation series we have tamed the 



■^This equation follows from noting that an infinitesimal change of variables 1 
leaves the integral J^ d(j)(j}e~^^'^^ invariant. 



strong- coupling behavior of perturbation theory. That is, all the higher order 
terms appearing in the expansion (0) are only C(l) at strong coupling. 

<S~So>=-\ + Oil/g) 
-l<{S-Sof>c,o=-^ + Oil/g) 



Indeed the first few terms of the reorganized perturbation series fall off nicely 
with the order. This provides an intrinsic reason to expect that truncating 
the reorganized perturbation series at low orders should provide a good ap- 
proximation to the exact result .0 



2.2 Varieties of gap equations 

Consider a set of degrees of freedom 0,, governed by an action S{(t>i). We 
wish to approximate this system with a simpler action SQ{(j)i). For the most 
part we will take 5*0 to be Gaussian in the fundamental fields. 



^^-J2^^ 



2^r 

But in some situations, in particular for gauge theories, we will be led to use 
a more general action, which is exactly soluble but not strictly Gaussian, of 
the form 



5o = 5:^a 



a 

Here the Aq are adjustable parameters, and the Oa are composite operators 
built out of the fundamental fields. 

In either case, the adjustable parameters appearing in 5*0 will be chosen 
to satisfy a set of gap equations. There are several prescriptions for writing 
down gap equations. All the gap equations given below have one feature in 
common: they correspond to critical points of an effective action.^ We now 
present several prescriptions for writing down gap equations. 

• Lowest-order Schwinger- Dyson gap equations 



''Although we haven't investigated the issue, it seems too much to hope that the reor- 
ganized perturbation series is actually convergent. 

^For gap equations which follow from a variational principle, the gap equations give 
a global minimum of the effective action. But in general the gap equations only give a 
critical point. 



For a generic degree of freedom 0j one has the Schwinger- Dyson equa- 
tion 

< 0.|^ >= 1 (9) 

d(pi 

(no sum on i). This is an exact relation among correlation functions 
in the full theory, which can be derived by the standard procedure 
sketched in section 2.1. It seems natural to choose the parameters 
of 5*0 in such a way that the same relation holds between correlation 
functions evaluated in Sq. That is, a natural set of gap equations 
corresponds to demanding that for each i 

This gap equation can be obtained from an effective action, as follows. 
Assuming that 5*0 is strictly Gaussian, and writing S* as a polynomial 
in the fundamental fields, one can show that the gap equations ([To|) are 
equivalent to requiring 

^ (/3Fo+<^-S'o>o)=0. (11) 



(9cr, 

That is, these gap equations make the first two terms of the perturba- 
tion series for (3F stationary with respect to variations of the af. 

Operator gap equations 

It is natural to ask that expectation values of the operators Oa agree 
when calculated in 5*0 and in S: 

< Oa > = < Oa >0 • (12) 

For example, when discussing gauge theories, we will introduce the 
Wilson loop operator O = TrPe*^^. Then the condition (|12D demands 



that Wilson loops agree when computed in 5*0 and in S*. As another 
example, suppose that Sq is Gaussian, with Oi = (pf. Then (|l^) means 
choosing the Gaussian widths af to be the exact 2-point functions of 
the full theory. 

As it stands (|12D is not a useful equation. But it can be rewritten as 

< ae-(^-^«) >c,o=< a >o (13) 

and expanded in powers of S* — 5*0. At leading order this condition is 
trivially satisfied. But at first order we get the gap equation^ 

< Oa{S- So) >c,o=0. (14) 



^It is tempting to try to improve this gap equation by expanding (|l3) to higher orders 
in 5 — 5*0. But this rarely seems to lead to useful gap equations. 



To relate this to an effective action, note the identities 
d 1 

-/3Fo = -^< a, >o 



9A " '' \l 



'a 



d I 1 

<S-So>0=^< Oa{S - So) >C,0 +T^ < Ca >0 • 



9A, ^^ A2 ^^"^" -u, ^.,u , ^2 



Thus the gap equation (jl^) is equivalent to requiring 

d 



dXa 



{f3Fo+<S-So>o)=0. (15) 



Note that if 5*0 is strictly Gaussian then these operator gap equations 
reduce to flTTf). 



Variational gap equations 

In some cases gap equations can be obtained from a variational prin- 
ciple. For purely bosonic theories one has a bound on the free energy, 
that 

PF < PFo+ <S-So>o ■ 

Thus one can formulate a variational principle, based on minimizing 



the right hand side of the inequality with respect to A^ [^ . Note that 
the resulting gap equations are equivalent to (|15|). 

The gap equations that we have discussed so far are all essentially equiv- 
alent: they follow from varying PFq+ < S — Sq >q. Unfortunately, they will 
not be sufficient to analyze many of the theories of interest in this paper. For 
example, suppose that 5*0 is purely Gaussian, and that S only has non-trivial 
3-point couplings. Such couplings do not contribute to /5Fo+ < S — Sq >o, so 
the gap equations discussed above are not sensitive to the interactions. But 
then, if some field is massless at tree level, ^Fq will suffer from an infrared 
divergence. 

To resolve these difficulties, we turn to the 

• One-loop gap equations 

The gap equation of section 2.1 resummed an infinite set of Feynman 
diagrams, involving one-loop self-energy corrections to the propagator. 
This can be extended to theories with both 3-point and 4-point inter- 
actions, as follows. Write S = Sjj + Sm + Sjy + ■ ■ ■ as a sum of 2-point, 
3-point, 4-point, . . . couplings, and assume that 5*0 is strictly Gaussian. 
Consider the quantity 

-^eflF = PFq+ < Sjj — Sq >o + < Sjy >0 " -^ < (Sjij) >c,0 ■ 



Q 



+ ^^^ + 



O 



Figure 2: One-loop gap equations. Heavy lines are dressed propagators and 
thin lines are bare propagators. 

Requiring -^les = gives a set of gap equations which resum all 
one-loop self-energy corrections to the propagators, built from both 
3-point and 4-point vertices. These gap equations are illustrated in 
Fig. 2. They are sufficient to cure the infrared divergences for most 
theories of interest. The quantity I^s can be identified with the 2PI 



effective action of Cornwall, Jackiw and Tomboulis [E^, as calculated 



in a two-loop approximation. 

In general it does not seem possible to give a prescription for choosing the 
best set of gap equations. Indeed it is not even clear what 'best' means - it 
may depend on what one wishes to calculate. For example, in some cases, the 
variational principle gives the optimal gap equations to use when estimating 
the free energy. A different set of gap equations will give a worse estimate 
of the free energy, but may give a better estimate of other quantities, such 
as 2-point correlations. The only general advice we can offer is that the 
gap equations must be chosen to cure all the infrared divergences which are 
present in the model. If possible, the gap equations should also respect all 
the symmetries of the model.0 



2.3 Difficulties with gauge symmetry 

There is a serious difficulty which one encounters when applying the Gaussian 
approximation to a gauge theory. The underlying gauge invariance implies 
Ward identities, which relate different Green's functions. But the Gaussian 
approximation often violates these identities, and this can lead to inconsis- 
tencies. 

To illustrate the problem, consider an action 5'(0j) which only has 3-point 
and 4-point couplings - for example ordinary Yang-Mills theory. We assume 
that the one-point functions < (pi > vanish. Then the exact Schwinger- Dyson 
equations < (pij^ >= 1 can be rewritten in terms of proper (IPI) vertices 



as shown in Fig. 3 |28|, ^, ^ 



''This becomes problematic when deahng with gauge theory, as we discuss in the next 







Figure 3: Schwinger- Dyson equations for a theory with 3-point and 4-point 
couphngs but no tadpoles. The sohd blobs are dressed propagators; the 
empty circles are IPI vertices. All external lines are amputated. 

The one-loop gap equations which we discussed in the previous subsec- 
tion are a truncation of this system. They correspond to dropping the two- 
loop contributions to the Schwinger- Dyson equations and approximating the 
dressed IPI 3-point vertex with a bare vertex. This truncation works well for 
many theories, as we will see. But in an abelian gauge theory, for example, 
a Ward identity relates the dressed 3-point vertex to a dressed propagator. 
This makes it inconsistent to work with dressed propagators but bare ver- 
tices; the inconsistency shows up in the fact that the one-loop gap equations 
with bare vertices do not have a solution.^ Methods for dealing with this dif- 
ficulty have been proposed [^, ^, but they do not seem to be compatible 
with manifest supersymmetry. An adequate analysis of the 0-brane quantum 
mechanics calls for a more satisfactory treatment of this issue [^. In section 
5, when we apply the Gaussian approximation to a bosonic large- iV gauge 
theory, we will dodge the issue of Ward identities, by avoiding the use of the 
full set of one-loop gap equations. 



3 Scalar models with supersymmetry 



In this section we apply the Gaussian approximation to study the strong- 
coupling behavior of some quantum mechanics problems having Af = 2 su- 
persymmetry but no gauge symmetry. The goal is twofold: to learn how 

subsection. 

^This result was obtained in conjunction with David Lowe. 



to apply the Gaussian approximation to a supersymmetric system, and to 
test the Gaussian approximation against some of the well understood strong- 
coupling dynamics present in these simple systems. The main lesson is that 
the Gaussian approximation works well, provided that it is formulated in a 
way which respects supersymmetry. Effectively this requires formulating the 
Gaussian approximation in superspace: gap equations must be introduced 
for the auxiliary fields. 



3.1 A model with a moduli space 

We begin by considering a model with three scalar superfields and a super- 
potential W = $i$2'^3- The Minkowski space action is (see Appendix A for 
supersymmetry conventions) 

Sm = — dt -^a<Pa + ilpalpa + -fafa " -faMcSabc + MblpcSabc 

where Sabc = +1 if «, b, c are distinct, otherwise. Non-perturbatively it is 
known that the model has a moduli space of vacua^ localized along the 0i, 
4>2, 03 coordinate axes. Thus the spectrum of the Hamiltonian is continuous 
from zero, and the free energy falls off like a power law at low temperatures. 

In studying this model, our goal is to show that the Gaussian approxima- 
tion captures this power-law behavior of the free energy at low temperatures. 
This is not a trivial result, because going to low temperatures is equiva- 
lent to going to strong coupling: the coupling constant g"^ is dimensionful, 
g"^ ~ (length)"'^, so the dimensionless effective coupling is g'^^ = g'^/T^- We 
henceforth adopt units in which g"^ = 1. 

We use an imaginary time formalism, Wick rotating according to 

Se = -iSM T = it faE = -ifaM ■ 

Note that we have to Wick rotate the auxiliary fields in order to make the 
quadratic terms in the action well-behaved in Euclidean space. Expanding 
the fields in Fourier modes 



Mr) = ^ J2 "^re 



i2-iTrT/f3 



Mr) = :i=E/r^ 



J27r/r//3 



^in the sense of the Born-Oppenheimer approximation 

10 



the Euclidean action is 

2 



pa 



+J73 E '>?A/nS.i« - 75 E '^f'/^vV.'s.fc. 

^ ^ l+m+n=0 ^ ^ i-r+s=0 

For a Gaussian action we take (the i?-parity symmetry discussed in Appendix 
A forbids any - / mixing) 

Next we need to choose a set of gap equations. This model only has cubic 
interactions, so a reasonable choice is the one-loop gap equations discussed 
in section 2.2. These read 



/27r/^ ^ 






?Tt+n=i r—s=l 



^2 = ^ + ^Y. -l-l (16) 

1 27rr 2 v^ 2, 

l+s=r 

where hr = —igr is real. Note that time-reversal invariance requires the 
symmetry properties 



a_, = a, T_, = T, 



Ifj-^r' I hr' 



Meanwhile the first few terms in the reorganized perturbation series for the 
free energy are 

I I r 



e((^)v-i)4ew-i)-3e(- 

3E((T'^-?-i)'-!E(-Nif 

3 / 2nr \^ 3 3 



— - < (S — So) >c,o= — 



l+m+n=0 /— r+s=0 
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Figure 4: I3F vs. (3 for the $i$2^3 model, as given by the sum of the first 
three terms of the reorganized perturbation series. Numerical calculations 
were performed at the indicated points. 
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Figure 5: Plots of (— ) individual terms in the free energy vs. (3 in the low 
temperature regime. The scale is log-log, and we shifted (3Fq by a constant 
to display its power-law behavior. The top curve is — < S — Sq >o, the 
middle curve is +| < {S — SqY >c,0) and the bottom curve is —{[3Fq + 1.04). 
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To solve the gap equations ([16|) one has to resort to numerical methods. 
An effective technique is to start by solving the gap equations analytically at 
high temperatures, where the loop sums are dominated by the bosonic zero 
modes crf^Q, rf^^. One can then repeatedly use the Newton- Raphson method 
[]35| to solve the system at a sequence of successively smaller temperatures, 



using the solution at inverse temperature (3 as the starting point for Newton- 
Raphson at /? + A/?. Further details on the numerical algorithm are given in 
appendix B. 

We have solved this system numerically up to /? = 28.5, corresponding to 
a dimensionless coupling gl^ = g^P^ ;^ 10^. The resulting free energy 

PF ^ (3Fo+ <S-So>o-l<{S- SoY >c,o (17) 

is shown in Fig. 4. From the plot it seems clear that the energy of the sys- 
tem - given by the slope of jSF - vanishes as the temperature goes to zero, 
so the Gaussian approximation has captured the fact that supersymmetry 
is unbroken in this model. Note that loop corrections to the auxiliary field 
propagators were crucial in obtaining this result - if we had formulated a 
Gaussian approximation just for the physical degrees of freedom, the ap- 
proximation itself would have explicitly broken supersymmetry. 

Moreover, in the low-temperature regime, the numerical results for the 
individual terms in the free energy are well fit by 

f3Fo = -1.04 - 0.705/3-^-^^ 

<S-So >o= -2.58(3-^-^^ (18) 

-^<iS-Sof>cfl=-l.77p-'-''. 

This behavior can be seen in Fig. 5, where we have plotted these quantities 
on a log-log scale. 

Evidently (3Fq approaches a non-zero constant at low temperatures, in- 
dicating an effective ground-state degeneracy of the system. It would be 
interesting to understand whether such a degeneracy is really present, or 
whether it is perhaps an artifact of the Gaussian approximation.^ 

What is more interesting for our purposes is that the subleading behavior 
of PFq, as well as the corrections < 5 — 5*0 >o and — ^ < {S — SqY >c,05 all 
vanish like a power law at low temperatures, with about the same exponent. 



^"The spectrum of the Hamiltonian is continuous from zero, so the usual argument that 
pF ^ as /3 — > oo does not apply. In fact it is not even clear whether the zero-temperature 
partition function for this model is well-defined, since it may depend on how one regulates 
divergences coming from the infinite volume of moduli space. 

14 



It seems plausible that all the higher-order corrections will also follow roughly 
the same power law. If this is the case then the Gaussian approximation 
provides a good estimate of the exponent. The coefficient in front of the 
power law is more difficult to obtain, but one might hope that keeping just 
the first three terms as in ([TtD is not such a bad approximation. This leads 
to a numerical fit 

(3F ^ -1.04 - 4.77/3-^-°^ . 



3.2 A model with a mass gap 

Next we consider a model containing a single scalar superfield with superpo- 
tential W = j^"^, corresponding to the Minkowski action 

Sm = \J dt ^02 + ,^^ + 1/2 _ ^^3 ^ 3^2^^ _ 

Non-perturbatively this theory is known to have a unique supersymmetric 
vacuum localized near = |]5B[ . In this section we study the finite tempera- 
ture behavior of this model using the Gaussian approximation, and show that 
it reproduces this known behavior in the limit of zero temperature. Moreover 
we will use the Gaussian approximation to obtain an estimate for the energy 
gap to the first excited state of the system. This is non-trivial because, just 
as in the previous section, going to low temperature is equivalent to going to 
strong coupling: the dimensionless effective coupling is gl^ = g^/T^- From 
now on we set g^ = 1. 

The Euclidean action in Fourier modes reads 

Q 

+ -n ^ (pl(t>m(pnfo - "^ ^ (t>l4>mAlps ■ 

^ l+m+n+o=0 '^ l+m-r+s=0 

For a Gaussian action we take 

Note that one has to allow the fields and / to mix in the Gaussian ac- 
tion, because the superpotential breaks the i?-parity symmetry discussed in 
appendix A. Thus we have the correlators 

< 4>l4>m >0 = O-/ Sl+m < 4>lfm >0= pA+m 

< flfm >0 = r^Sl+rn < ^r^s >0= -fi'r'^rs ■ 

15 



For this model the one-loop gap equations are identical to the lowest-order 
Schwinger-Dyson gap equations. In either case one obtains the system 



PI rf) - V i; + v fE™^^ 

1 2-nr 3 v-^ 9 / x 

The solution to the gap equations has the form 



O", 



(2nl/PY + ml 



Pi = -irrifaf 

2 1 2 2 

Ti = l-m^ai 



9r 



—i2nr/l3 — nif 

parameterized by effective bose and fermi masses tti^, rrif. It is convenient 
to introduce the size of the state 

^rms = <(t>{rf> 



^E-'^ (20) 



^ I 



2mi, tanh(/5m;,/2) 

The gap equations then imply consistency conditions which fix m^, rrif as 
functions of the temperature. 

ml = 27R^^^-3tanh{3f3R^^„j2) 

mf = 3Rl,, (21) 

Meanwhile the first two terms in the expression for the free energy are 
PFo = -^ E log {aM -pD+IY: log \9r\' 

I r 

< S - So >o= i E ((^)V? - i) + 1 ^(.f _ 1) + ^ (,!^<,., + 1 

I m I r 

With the help of the gap equations this reduces to 

(3Fq + <S-Sq>q = log(2 smh{pmb/2)) - log(2 cosh(/5m//2)) 

1 9 

-^f^^bRrms + 2/^-^rms • (22) 

16 



The nice thing is that if you expand for low temperatures (equivalently 
strong couphng - recall that the dimensionless coupling is g'^f^ = g'^/T'^) then 
the solution to the consistency conditions (pO|), ( pT]) is 



Prrih 



y|/3 + v^/?e-^v^ + 0(e 



-fsVe 



The corresponding free energy (E^) is 



(3Fo + <S-So >o= -2e~'^V^ + O (e"^^) . (23) 

But this is exactly the low temperature expansion one would expect for 
the free energy of a system with a unique zero energy ground state (a 'short 
multiplet' of the Af = 2 supersymmetry) plus a degenerate pair of first excited 
states with energy a/3/2 (a 'long multiplet' of A/" = 2). 

We conclude by examining the second-order corrections to the free energy 
of this model. After imposing the gap equations one finds that 



1 / n n \2 ^ n "^ ^ \~~^ 2 2 i-^iTTl-c ^-— v 2 2 2 2 

'^ Z+m-r+s=0 '^ Z+m+ra+o=0 



The three-loop sums are best evaluated in position space, using (0 < r < /5) 

1 V^ 1 -i2-.iTip ^ coshmfc(r-/J/2) 

P ^ (27r///3)2 + ml 2mb sinh(/3mfe/2) 

1 ,r^ 1 , p"1/(t-/3/2) 



E^ — » 



13 ^ -i2Tir/l3 -ruf 2 cosh(/5m//2) ' 

to obtain 

1 2. 1^^3 9{/3mh + smh/3mk) 

< (O — On) >C0= —plfTif ;7^ 

2 ^ "^ ^'" 9^ f 16m3(l + cosh/3m/)sinh2^ 
3/5"^f ,3^ 1 . , ^ 1 



V 



/5mfe H — sinh [3mh -\ sinh 2(3mh) . 



Amlsmh^^'s"^ 2 "^ " 16 

At low temperatures this reduces to 

-l<iS- So? >c,o= -2y|/3e-^v^ . (24) 

Although exponentially suppressed, this larger than the contribution of the 
first two terms ( ^3]) by a factor of /3. Combining the results (^3]), ( ^4]) we see 
that the expression for the free energy at low temperatures has the form 

-I3V6 



pF = -26-^^/^ (l + 13^/2/3 + ■■■') +o(e 
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Presumably the second-order term can be interpreted as correcting the energy 
gap to the first excited state. It would be interesting to find a set of gap 
equations for which the second order term in the free energy is subdominant 
at low temperatures. 



3.3 A model which breaks super symmetry 

Finally, we consider a model with a single scalar superfield and a superpo- 
tential W = |$^, corresponding to the Minkowski action 

g^ J 2 2 

The dimensionless coupling constant is g^f^ = g'^/T^, and we henceforth set 
9' = I- 

This model is known to break supersymmetry [Q. This is an effect which 
cannot happen at any finite order in perturbation theory, but we will see that 
it is captured by the Gaussian approximation: the ground state energy comes 
out positive. This provides an interesting contrast to the usual analysis of 
supersymmetry breaking via instantons [0, which can be applied to this 
model once a sufficiently large mass perturbation is introduced. A related 
analysis of supersymmetry breaking in the 0{N) cr-model was carried out by 
Zanon [^ . 

The 0(2) R-symmetry prohibits any - / mixing, and implies < >= 0, 
but allows a tadpole for /. The expectation value of / is the order parameter 
for supersymmetry breaking. To take it into account we shift / — >■ a + /, 
with < / >= 0. The Euclidean action i^ 



5- = lECfX^'^-'+lEf'f-'-i: 

I ^ ^ '' I r 

--/?a^ - ia\/^fo + a ^ 0;0_i 



i—rAA 



+ -7^ Y^ fl(pm(pn - ^ XI ^I'^r'ips ■ (25) 

'^ ^ l+m+n=0 ^ ^ l-r+s=0 



For a Gaussian action we take 






^We Wick rotate the fluctuation / but not the expectation value a. 
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We will use the one-loop gap equations, which read {hr = —igr 



' \ ' / ' fn+n=l ' r-s=l 

:^ = 1 + 1 E -^-n (27) 

l+s=r 

The value of a is fixed by requiring < / >= 0. At the one loop level this 
gives 

« = ^E^'- (28) 

Note that in this approximation the order parameter for supersymmetry 
breaking is the spread of the ground state wavefunction! We approximate 
the free energy as I3F ^ I3Fq+ < S — Sq >o, with 



I I r 

<s-s„>„4i:((^)w-i)4ew-i)-i:(^^-i) 



The gap equations for this model must be solved numerically. We show 
the result for a =< / > in Fig. 6. The key point is that < / > approaches a 
non-zero constant at low temperatures, < / >~ 0.50, so that supersymmetry 
is indeed broken. This is refiected in a positive slope of f3F, corresponding to 
positive vacuum energy. The numerical results can be fit at low temperatures 
by 

pFo = -0.33 + 0.21/3 
<S-So >o= -0.13 + 0.057/5 

pF^ -0.46 + 0.27/3. 

Note that, although the vacuum energy is non-zero, it is not given by the 
classical formula ii^ = |</>^-we are at strong coupling and quantum 
corrections are important. 

To conclude, let us note that it is not really necessary to solve the full 
set of one-loop gap equations to study this model. Rather, it is a reasonable 
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Figure 6: Auxiliary field expectation value vs. f3 in the $^ model. 
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approximation to truncate the gap equations (127|) to 



1 


/27iiy 


^f 


\ (3 ) 


1 


= 1 


1 


2nr 


hr 






+ 2a 



while keeping equation (pS]) for a. In this approximation the gap equations 
are solved by 

1 



a 



8tanh'^{(3^/a/2) 
while the free energy is given by 

f3F ^ (3Fo+ <S-So >o= logsinh(/?v^) - -/?«^ • 
Note that at low temperatures a -^ 1/2 and PF -^ |/5. 



Gauge dynamics and the Gross- Witten tran- 
sition 



In this section we discuss a generic feature of large-A^ gauge theory: the 
existence of a Gross- Witten transition [^.Q We will argue that such a 
transition separates the weak-coupling and strong-coupling regimes of the 
gauge theory. 

Suppose we have a -|- 1 dimensional U{N) gauge theory at finite tem- 
perature. Although the gauge field Aq{t) has no local degrees of freedom, 
one can construct a non-local observable by diagonalizing the Wilson line. 

/ e*"i \ 

y e*"^ y 

The eigenvalues ai, . . . , aN characterize the gauge- invariant degrees of free- 
dom contained in Aq (up to permutation by the Weyl group). 



As discussed in |^, for most theories the eigenvalue distribution will 



be quite different at high and low temperatures. This can be understood 

^^By Gross- Witten transition we mean a large- A^ phase transition triggered by a change 
in the eigenvalue distribution, not necessarily a third order phase transition. 
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as a competition between entropy and energy. The Haar measure for the 
eigenvalues, corresponding to Wilson lines which are uniformly distributed 
over the U{N) group manifold, is given by [|^ 



Y\ dai n (1 - cos{ai - aj)) 



i i<j 

The measure vanishes whenever two eigenvalues coincide. So in the absence 
of any other interactions, the eigenvalues tend to repel each other, and spread 
out over the circle in order to maximize their entropy. 

On the other hand, energetic considerations often lead to attractive in- 
teractions between eigenvalues. As a simple example, suppose we minimally 
couple the gauge field to a massive adjoint scalar field. Integrating out the 
scalar field produces an additional contribution to the measure for the eigen- 
values, given by 



l[det ' 



i<j 



idr + -{ai-aj)] +m^ 



n 



. 2 (cosh/5m — cos(Q;i — a^)) 



(29) 



In the high temperature regime /3m <^ 1 this gives a strong attraction be- 
tween eigenvalues, and the eigenvalues will tend to cluster to minimize their 
energy. But at low temperatures the potential p^ ) flattens out, and entropy 
wins: the eigenvalues will spread out around the circle. 

This phenomenon seems to be quite generic: eigenvalues will tend to clus- 
ter at high temperatures, and spread out at low temperatures. This should 
happen even in supersymmetric theories.|3 Of course in quantum mechanics 
at finite N one has smooth crossover between these two qualitatively different 
behaviors. But in the large A^ limit a sharp Gross- Witten phase transition 
should occur when the eigenvalues can first explore the entire circle. 

We conclude with a few remarks on this phenomenon in the context of 
D-brane physics; for a review of the phases of D-brane gauge theories see 



0. 



1. The phenomenon is well-known, at least for D-branes which are wrapped 
on a spatial circle. It is responsible for the transition between singly- 
wound and multiply- wound strings discussed in |H2|, H3| , and the related 



black hole transitions discussed in [44|. We've merely pointed out that 



^^While bosons generate an attractive potential which is minimized when ai — aj — 0, 
fermions give rise to a repulsive potential that repels away from ai — aj = n. 
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at finite temperature the same transition takes place in Euclidean time: 
0-brane worldlines can become multiply wound around the time direc- 
tion. Only the multiply-wound phase corresponds to a black hole. 

2. In a T-dual picture the eigenvalues of the holonomy become the po- 
sitions of D-instantons on the dual circle. At high temperatures the 
dual circle is large, and the D-instantons tend to cluster. But at low 
temperatures, when the dual circle is small, the D-instantons will de- 
localize. Note that only in the low-temperature regime could one hope 
to describe the configuration as a black hole, using a static solution to 
the supergravity equations. 

Let us make a crude estimate of the transition temperature from the 
D-instanton point of view. Consider a system of A^ D-instantons in 
IIB string theory with string coupling qb and string length /<,. Rea- 



soning analogous to [^ ^ shows that in uncompactified space the 
D-instantons will cluster and fill out a region of size L ~ [gBNY^'^ls. 
Compactifying the time direction on a circle of size (3bi we expect the 
D-instantons to delocalize if j3b < L. Mapping to type IIA via 

(3a = a' /(3b g\ = gla'/(3l 

we conclude that the 0-brane worldlines are multiply-wound when 

(3a > (3crit = {9aN)-'/% = {qImN)-^'^ . (30) 

The transition at /5crit can be seen in several ways from the supergravity 
point of view. The background geometry of a system of 0-branes in the 
decoupling limit is well-defined above (3cr\t- Below /3crit the curvature at 
the horizon exceeds the string scale, and the dual gauge theory becomes 
weakly coupled [p4| . 

The transition can also be identified with the Horowitz-Polchinski cor- 



respondence point [0. Consider a non-extremal black hole in ten di- 



mensions with 0-brane charge. Following |T^, WM , one can map this to 



a ten-dimensional Schwarzschild black hole by lifting to eleven dimen- 
sions, boosting in the covering space, and recompactifying. We discuss 
this transformation in Appendix C. If the original charged black hole is 
at the critical temperature, then the equivalent neutral black hole turns 
out to have a Schwarzschild radius close to the string scale. That is, 
the temperature Tcrit maps to the Horowitz-Polchinski correspondence 
point, at which a black hole is on the verge of becoming an elementary 
string state. 

4. Although our primary interest is in quantum mechanics, it is natural 
to suppose that similar transitions happen in other dimensions. For 
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example, consider 3 + 1 dimensional A/" = 4 large-A^ Yang-Mills on 
M^ X 5*^. For any non-zero temperature it seems natural to expect that 
a Gross- Witten transition occurs when the 't Hooft coupling A = Qym^ 
is of order one. 

To conclude, we expect that Gross- Witten transitions are generically 
present in large-A^ gauge theories at finite temperature. In theories with a 
dual string interpretation this suggests that a phase transition occurs, which 
separates the non-geometric phase from the phase with a well-defined back- 
ground geometry (see also [^, ^). 



5 Large-A^ models in the Gaussian approxi- 
mation 



In this section we use the Gaussian approximation to study some large-iV 
matrix quantum mechanics problems. One motivation for the problems we 
consider is that they arise as subsectors of the 0-brane quantum mechanics. 
The 0-brane quantum mechanics can be written in terms of A/" = 2 superfields 
as |5^ (see Appendix A) 

9ym J ^ ^ ^ "J ■' 

(31) 

Here JFj is the field strength constructed from a gauge multiplet Fq,, the $a 
are a collection of seven adjoint scalar multiplets transforming in the 7 of 
G2, and fabc is a suitably normalized^ totally antisymmetric G2-invariant 
tensor. 13 

In the next section we drop all the fermions, and (with nine scalar fields) 
study the bosonic sector of this action. In the following section we drop the 
gauge multiplet Fq,, and study the supersymmetric dynamics of the scalar 
multiplets $„• 



Jabcjabd — 2 cd 



"^^Under an SO{2) x 5*0(7) subgroup of the 5*0(9) i?-syminetry the physical fermions 
transform as an 8 of 50(7) , which decomposes into 7 ® 1 of G2 . Thus one of the fermions 
becomes part of the M — 2 gauge multiplet, while the remaining seven transform in the 7 

0fG2. 
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5.1 Bosonic large- A^ gauge theory 

In this section we study a purely bosonic quantum mechanics problem: the 
reduction of rf + 1 dimensional U{N) Yang- Mills theory to + 1 dimensions. 
We will study this theory at leading order for large A^ using the Gaussian 
approximation. The goal is to show that the Gaussian approximation re- 
spects 't Hooft large-iV counting, and can capture the Gross- Witten phase 
transition which we expect to be present. We will also learn a few useful 
lessons about 0-brane quantum mechanics. 

The degrees of freedom of this model consist of d adjoint scalar fields Xj, 
coupled to a gauge field Aq. The Euclidean action is 

Se = ^ [ drTrllD^X'D^X' -hx\X^][X\X^]\ (32) 

9ym ^ I ^ 4 J 

where D^ = d-r + i[-Ao, ■]. We set drAo = 0, so that Aq = const. = -^Aqq. 
Corresponding to this gauge choice we must introduce a system of periodic 
ghost fields a, a. In Fourier modes the complete action is0 

^ E ^Tr(X;[Aoo, XU]) + ^1^ Yl ^Tr(a,[Aoo, a,]) 



^^MV^Y ^ ^^MV^I^ P 



-^^Tr([Aoo,X;][Aoo,Xl,]) - -^ Y. Tr([X;,X4][X;,X^]) 
-^UyaiP , ^9ymP 



l+m+n+o=0 

Next we need to choose an action 5*0 to use as the basis for the approxima- 
tion.|^ Up to now we have always taken 5*0 to be quadratic in the fundamental 
fields. But this is not appropriate for gauge theory, because the eigenvalues 
of Aqq are angular variables. So it is better to work with the holonomy 

The 'Gaussian' action we propose to use is 

S, = -jTiiU + U^) + Y ^Tr(X;XlJ - J^ jMc^i^^i) ■ 

Strictly speaking, we propose to use this action to describe the relative eigen- 
values of the U{N) gauge theory. That is, U really takes values in the group 

^^Note that our gauge condition leaves a residual global U{N) gauge symmetry unbro- 
ken. We have suppressed the corresponding ghost zero mode a;=o- 

^^Out of habit we will keep calling it the Gaussian approximation, even when 5*0 is not 
strictly Gaussian. 
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SU{N)/Zn. But we will largely ignore this distinction, since it does not 
affect leading large- A^ counting. 

Although the action 5*0 is not strictly Gaussian, it is soluble, at least in 
the large- A^ limit: the action for U is the one-plaquette model studied by 



Gross and Witten [^. We briefly recall some results. The one-plaquette 
partition function is 

\ -A^VA^ A > 2 

with a third-order phase transition at A = 2. Diagonalizing the holonomy, 
U = diag.(e*°\ . . . , e*°^), one has the distribution of eigenvalues 



2 „^„ a I \ „;^2 

pn[a 



.Acosfjf-sin^f A<2 



^(1 + fcosa) A>2 

The expectation value of the Wilson loop is 

JN{l-X/4) A<2 

< TtU >□= < 

[ N/X A > 2 

Finally, the 2-point correlators in 5*0 are given by 

< ^00 AB^OQ CD >0 = Po^AD^BC 

< ^lAB^LcD >0 = (rf6'^6i+mSADSBC 
< OilABOimCD >0 = Sidim^AD^BC 

where 

^[li^(l-|) + (l-!)log(l-f)-l] A<2 

A> 2 






involves a dilogarithm. 



The first two terms in the expansion of the free energy are given by 
2 



PFo = PFni\)-^Yl l°g ^' + ^' E l°g ^' (33) 



- ni \ nf A 
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In writing these expressions we have kept only the leading large- A^ behavior, 
coming from planar diagrams. 

The next step is to choose a set of gap equations. For this system we 
adopt the operator gap equations discussed in section 2.2. That is, we will 
require that PFo+ < 5 — 5*0 >o be stationary with respect to variations of 
the parameters A, af, si. 

Let us first consider the gap equation from varying af, which reads 



^? 9ym V P 



) ^91mP''^ 91mP V ' 



1 _ 1 /27r/Y ^ 2iV ^2 , 2iV(rf 



The solution is af = 5'y^^/((27r///9)^ + rn?), parameterized by an effective 

thermal mass 

2N 
m' = —pl + 2{d-l)Rl^,. (34) 

Here we have introduced the 'size' of the system, defined by the range of 
eigenvalues of the matrices Xi . 

-^rms = TrTr < {X\t)Y > (no sum on i) 

91mN (35^ 



2mtanh(/3m/2) 
Next we consider the gap equation for the ghost propagator, which reads 

1 1 f2TllV , , 

si 9ym \ P J 

Thus in this approximation the ghosts are free and decoupled. At leading 
order they make a trivial contribution to the partition function, since in + 1 
dimensions one has det' 9,- = 1 . 

Finally, we have the gap equation from varying A, which reads 

1 - (1 - f ) log (i^) = gl^,Np/AdR^,^ A < 2 

[6b) 
X = g'yMN(3/2dR^^, A>2. 



By making use of the gap equations, the expressions (|33| ) for the free 
energy can be simplified to 

f3Fo = f3Fu{X) + A^^t/log(2sinh(/3m/2)) 

< g - g„ >„= g^ < TM >. - ^^5' " ^' (RL.f . (37) 
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Before studying these equations further, let us note that this approxima- 
tion respects 't Hooft large- iV counting. All the factors of Qym ^^^ ^ which 
appear in these equations are dictated by 't Hooft counting plus dimensional 
analysis, and they can be eliminated by introducing the following rescaled 
dimensionless quantities: 

A = A 

pI = NpIMmNY/' af = Naf/{g'y,,Ny/' 5, = iVs./(^^,,iV)V3 

^ I 

J3F = ISF/N"^ . (38) 

The 't Hooft coupling Qym-^ ^"^^^ as a unit of (energy)'^. The effective di- 
mensionless coupling is then gyf^jN/T^, so the system is strongly coupled at 
low temperatures. 



The system of equations (0), (p5D, ( pq) suffices to determine A as a 
function of the temperature. The results for the case of rf = 9 scalar fields 
are plotted in Fig. 7. Note that a Gross- Witten transition occurs at /? = 9.0, 
the temperature at which A passes through 2. The approximation suggests 
that the transition is weakly second-order: the second derivative of (3Fq+ < 
S — So >o jumps by -0.0053 across the transition.[^ But this may be an 
artifact of the approximation; a third order phase transition would seem 
more plausible. 

At high temperatures (T > [gyM^Y^'^-, ^ i^i dimensionless terms T > 1) 
the Gaussian approximation gives the behavior (with /3F ^ /?-fo+ < 5' — 
'S'o >o) 



A 



1/2 



\d{2d-l) ' ^ 



(3F ^ ^N\d -1) log p (39) 

E ^ ^N\d-1)T 

This is the expected result in the high temperature regime: the N'^{d — 1) 
physical degrees of freedom making up the matrices X* are excited, with an 
energy per degree of freedom of order the temperature. 

^^The fact that the first derivative is continuous across the transition is guaranteed by 
the form of the gap equations that we are using. 
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Figure 7: A vs. (3 for bosonic gauge theory with 9 scalar fields. The temper- 
ature is measured in units of the 't Hooft coupling. 



At low temperatures (T < (gyM^Y^^^ or in dimensionless terms T < 1] 
the Gaussian approximation gives 



it. 



d 

2(d- 1)1/3 
3 



f5F ^ -N^d{d-lf'^p 
8 

3 



(40) 



E ^ lN'd{d~lY'\g'y^,N) 



1/3^^2 __ArN 1/3 



The physics here is equally simple: the N'^d oscillators have all been frozen 
into their ground state. We see that the ground state energy of an individual 
oscillator is very large, of order ((yfy^A^)i/^, while the ground state wave- 
function of an individual oscillator is extremely narrow, being a Gaussian of 
width 



It is only when all of the fiuctuations are added together in R^ras that the 
system as a whole has a very large ground state size, i?rms ~ (S'yA/^)^^^- 



Before arguing that the Gaussian approximation really does give reliable 
results for this system, let us pause to compare the behavior of this bosonic 
system to the behavior expected of 0-brane quantum mechanics ESl 061 . The 
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full 0-brane system corresponds to d = 9 and differs by including sixteen ad- 
joint fermions. In the high temperature regime, we expect that the fermions 
can be neglected on general grounds (from the path integral point of view 
they are squeezed out by their antiperiodic boundary conditions). So the 



results ( pOD should be a good approximation to the behavior of 0-branes at 



high temperature. On the other hand, in the low temperature regime the 



0-brane system has a dual supergravity description p4| . The size of the re- 
gion in which supergravity is valid ~ (gyM^Y^^ is expected to match the 
size of the ground state wavefunction of the supersymmetric gauge theory. 
We see that the scale {gyM^Y^^ already appears as the size of the bosonic 
ground state wavefunction. So it seems that including the fermions should 
not significantly change the size of the wavefunction. The fermions do en- 
ter in a crucial way in the supersymmetric theory, however, because they 
are necessary for cancelling the enormous bosonic ground state energy which 
appears in 



We will now argue that the Gaussian approximation provides a good de- 
scription of this bosonic gauge theory for all values of the temperature. Of 
course the Gaussian approximation is good in the high-temperature regime, 
where the gauge theory is weakly coupled and the action (^2D is approxi- 



mately Gaussian. But it is also a good approximation at low temperatures. 
To establish this, let us examine the second-order corrections to the ground 
state energy of the system. 

The second-order correction — ^ < {S—SqY >c,o is a sum of several terms. 
But contributions involving the gauge field propagator can be neglected at 
low temperatures, because pi —^ as (3 ^ oo (the zero mode of the gauge 
field becomes frozen as the circle decompactifies). Most of the remaining 
terms cancel when the gap equations are imposed, leaving just 

— - < ("S* — 5*0) >c,o 



2 



N^d{d-1) ^ 2222 

g^j^N^d{d-l)p /3^ 1 . w^ N 1 • wo^ . 

■^^^ ^^^ — j^ ( -f3m + - smh(/3m) + — smh{2(3m] 



16m5 sinh^(/?m/2) \8' 2 '' ' 16 

(the three- loop integral is easiest to do in position space). Combining this 
with previous results gives the first three terms in the reorganized perturba- 
tion series for the free energy, which in the low temperatures regime readsj^ 

1 1 



pF = -N^d{d - ly/'^p 



1- 



4 16{d-l] 



^^Presumably this model admits an expansion in 1/d. It would be interesting to under- 
stand this expansion, and its relation to the Gaussian approximation. 
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As claimed, the higher order terms in reorganized perturbation theory make 
a small correction to the free energy, even in the strong-coupling, low tem- 
perature regime. Based on the size of the second order correction, we can 
claim to have computed the ground state energy of this system to within 
roughly one percent accuracy for the case of 9 scalar fields. 



5.2 A supersymmetric large- A^ scalar model 

Consider a set of A^ x A^ Hermitian scalar multiplets $a, governed by an 
action which is a truncation of the 0-brane quantum mechanics (0). 

^= ^ /"rftd2^Tr|-^D-$,Z}„<l>,-^/„,,$J<l>5,<l>e]} (41) 

Here a = 1, . . . , 7 is an index in the 7 of G2 and fabc is a G2-iiivariant tensor. 
We will subsequently set Qym = 1- Fo^ ^ Gaussian action we take 

^0 = E ^Tr <Ptr-i + E ^Tr /,"/!, + $^ ^Tr ^."^ . (42) 

The one- loop gap equations are evaluated in the large A^ limit, by keeping 
only the planar contributions. We find that the gap equations have a familiar 
form 

' '^ '^ m+n=l ' r-s=l 

1 A 



P 



1 + - V A'A^ 



m+n=l 



1 27rr 2A v^ ,„ 

l+s=r 

where h^ = —igr is real and A = 14N. 

In fact, in the Gaussian approximation, this G'2-invariant model is iso- 
morphic to the $i$2$3 model studied in section 3.1. The gap equations (|i3| ) 
can be exactly mapped to those of the $i$2'^3 model by defining 



Pnew 


= A^^'P 


A,2 

-'new 


= A'/'A' 


e2 = 

^new 


= e' 


'''now 


= A'/'h. 



Restoring units, note that A is essentially the 't Hooft coupling, so that jS^ew 
is the physical temperature measured in units of the 't Hooft coupling. With 
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this redefinition one finds that PFq, < S — Sq >o and < {S — SqY >c,o for the 
G2 model transform into the corresponding quantities for the $i<l>2*^3 model, 
up to an overall factor of 7N'^/3. Thus all the results of section 3.1 can be 
taken over and applied to the G2 model. In particular, the free energy falls 
off as (3~J^^ at low temperatures. 



6 Conclusions 



In this paper we've discussed an approximation scheme which can be used to 
study strongly-coupled quantum mechanics problems. We've shown that the 
approximation can be formulated in a way which respects supersymmetry, 
by introducing gap equations for the auxiliary fields. We've also shown that 
it can be applied to large- iV theories: indeed 't Hooft scaling is automatic, 
provided that one keeps just the planar contribution to the gap equations. 

We used the approximation to study a number of toy models of scalar 
supersymmetric quantum mechanics, and found that it captures the correct 
qualitative behavior of the free energy at strong coupling. Moreover, it pro- 
vides a quantitative estimate of certain quantities, such as the ground state 
energy or low-lying density of states, which would be difficult to obtain by 
any other means. 

The approximation has difficulty dealing with gauge symmetry, since the 
gap equations tend to violate the Ward identities. Nonetheless, we were able 
to study a large-A^ gauge theory: the bosonic sector of 0-brane quantum 
mechanics. The Gaussian approximation provides a remarkably accurate de- 
scription of the ground state of this theory. The approximation also predicts 
that a Gross- Witten phase transition occurs when the effective coupling is of 
order one. 

We argued that such a phase transition should be present in a generic 
large-iV gauge theory, in particular for the full 0-brane quantum mechan- 
ics. This means that the perturbative gauge theory regime of the 0-brane 
quantum mechanics is separated from the supergravity regime by a mild 
phase transition. It also suggests that the Horowitz-Polchinski correspon- 
dence point is associated with a phase transition. 

Of course our motivation for this project was to develop techniques for 
understanding 0-brane quantum mechanics. We have been able to incor- 
porate many of the necessary ingredients into the approximation, including 
supersymmetry and large- iV counting. But an adequate analysis of the full 
0-brane quantum mechanics requires a more satisfactory treatment of gauge 
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invariance ||3^ . 
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A Supersymmetry in quantum mechanics 

In this appendix we review the superspace construction of A/" = 2 theories in 
+ 1 dimensions. 

With Af = 2 supersymmetry we have an 5*0(2) R-symmetry, with spinor 
indices a, j3 = 1,2 and vector indices i,j = 1,2. The S0{2)ji Dirac matrices 
are real, symmetric, and traceless; for example one can choose 7^ = ui and 
7^ = (73 where ui and a 3 are Pauli matrices. Given two spinors tpa and Xa, 
besides the invariant ipaXa, one can construct a second invariant which we 
will denote 

Af = 2 superspace has coordinates (t, 9a) where 9a is real.0 Denoting 9'^ = 
9°'9a, we normalize J (P9 9'^ = 1. 

We introduce the supercharges and supercovariant derivatives 



which obey the algebra 



Qa = 


d 

d9a 


-< 


Da = 


d 
d9a 


--4 


{Qa,Q(3} 


^ 


2i5ai3dt 


{Qa,Dp} 


= 





{Da, Dp} 


= 


-2i6apdt . 



^"Our convention is that complex conjugation reverses the order of Grassmann variables. 
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The simplest representation of supersymmetry is a real scalar superfield 
The general action for a collection of scalar superfields $a is 



r J 2 



^ = — 

9 

g^ J 2 2 (90a <90a<90fe 

where in the second line we have taken complex combinations ipa = -y^ (V^ia + 'i'^2a), 

In some cases we can extend the 5*0(2) /^-symmetry to an 0{2)ji, where 
the extra Z2 ('i?-parity') acts as 9^ -^ lap^^p- Both the measure d^O and the 
kinetic term are odd under i?-parity. Thus i?-parity is a symmetry if the 
superpotential is odd in $, with $ -^ — $ under R-parity. 

To describe gauge theories we introduce a real connection on superspace 

V„ = D^ + V^ 
Vt = dt + iVt 

with the gauge transformation d^Ta = «Vq,A, 6ATt = VfA. The correspond- 
ing field strengths J-'a/s, ^a are real superfields given by 

{Vi,V«} = To, 

We impose the so-called conventional constraint dapTafi = 0, which implies 
that Tap is a vector of S0{2)r. 

Tap = '^lapTi . 

This fixes Tt in terms of the spinor connection Fq,, but does not constrain 
Ta, as the only independent Bianchi identity 

V.^, = -tYap^p (44) 

is satisfied for any Ta- 

We write the expansion of Fq, in 'linear' components as 

The fields X* are physical scalars, while Aq, are their superpartners, d is an 
auxiliary boson, Xa are auxiliary fermions, and Aq is the 0+1 dimensional 
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gauge field. These linear component fields are of primary interest: by using 
them to formulate the Gaussian approximation we are effectively making a 
Gaussian approximation in superspace for the superfield Fq,, and such an 
approximation is guaranteed to respect supersymmetry. But it is also useful 
to define the 'covariant' component fields 

V* I — T« I 

-^ |cov = -J \0=Q 

1 • 

\ I ^ ^ I 

^q|cov — •'a 1 9=0 

\ t % • ■ % % 

= A„ + -Xa + -[^0, Xa] - ^lap[X\ Xp] + ^^afild, Xfi] - ^[Xa, XpXp_ 

^olcov = rt|e=o 
1 



— ^0 + ^XaXallin- 

These covariant component fields transform in the usual way under conven- 
tional gauge transformations, and are neutral under gauge transformations 
which depend on 9. Written in terms of covariant component fields the SYM 
action takes a simple form, 

SsYM = J- /dtrf2^Tr|-iv"^,V„^, 



9ym 



i 1 . . ^ 

+-KDoK - ^KYa(3[^\ A/3] I 



where Dq = dt + iAo. In terms of the linear component fields the action 
would only take this form in Wess-Zumino gauge. We must also introduce a 
ghost multiplet, for which we adopt the component expansion 

C = a + pJa + -f9^ 

where a and 7 are complex Grassmann fields and f3a is a complex boson. To 
the SYM action we add the ghost and gauge fixing terms 






"Jg.f. ~r »Jghost 2 

A'ym 
Note that the gauge-fixed action respects i?-parity. 



B Solving gap equations numerically 

The gap equations are an infinite set of coupled non-linear equations which 
determine the Fourier modes of the propagators. In general we must resort 
to numerical methods in order to solve them. 
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An important simplification follows from the fact that quantum mechanics 
is UV free. Thus at large momenta the propagators have free-field behavior, 
and at any given value of the temperature we only need to solve the gap 
equations to determine a finite number of Fourier modes. In the moduli 
space model of section 3.1 we obtained good results by solving for modes 
with -1.8(3 <l,r < 1.8/3. 

We are left with a finite set of equations, which we solved numerically 
using the Newton-Raphson method. Let us briefly summarize the method. 
Consider a system of equations 

Fi{xi,...,Xn) = i = l,...,n. (45) 

We start by guessing a solution XqW, and expand F to first order around Xom- 
Solving the linearized equations gives an improved guess 

^new = ^oid + 5x (46) 

where 

(47) 



ox — J -TqI^j Jjj 



dxj 



aiold 

If the initial guess is sufficiently close to a root, this procedure can be iterated 
a few times and will rapidly converge. We found that five or six iterations 
gave excellent results in the moduli space model. The Newton-Raphson 
method does suffer from poor global convergence, but this can be improved 



using a procedure known as backtracking |35 . 



Even with backtracking one needs an initial guess that is sufficiently close 
to the actual solution. One strategy is to start out at high temperatures. 
Then the gap equations are dominated by the zero modes, and one can obtain 
an approximate solution analytically. For example, consider the moduli space 
model. At high temperatures (/5 ^ 1) the gap equations are approximately 
solved by 






- v(3 at^, = {P/2niy 
tI = 1/2 r,% = l (48) 

hr = (3/27rr 

This can be used as a starting point for Newton-Raphson. Then one can 
solve the system at a sequence of successively lower temperatures, using the 
solution at one temperature as the initial guess for the next temperature. In 
the moduli space model we found that increasing /? by a factor of 1.2 at each 
step worked well. 

The gap equations can also be solved to determine the large-momentum 
behavior of the propagators at any temperature, up to constants that can be 
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extracted from the numerical solution for the low-momentum modes. For ex- 
ample, in the moduli space model, the propagators at large momenta behave 
as 

1 



a 



2 



T, 



' (27r///3)2 + m? 

' (27r///5)2 + m2 

27rr//3 



where 



,2 



2 °° 



^ l=-oo 

Note that the asymptotic forms of the propagators are related in the expected 
way, given by the naive supersymmetry Ward identities. 

In fact the high momentum behavior is important for computing the free 
energy. Even in quantum mechanics one has to deal with the divergent 
quantity 'YlA'^^^f- This sum, which appears in /3Fo, can be defined by 

oo 

^ log [af ((27r///?)2 + m^)] - 2 log(2 sinh(/?m/2)) . (50) 

— oo 

The modified sum converges, so it can be calculated numerically. 

C Black hole — SYM correspondence 

We discuss the connection between ten dimensional Schwarzschild black holes 
and finite temperature SYM quantum mechanics. We relate the two theo- 
ries by boosting an eleven dimension black string and recompactifying it 



on a smaller circle to keep the entropy the same [|12| , \l% . Although this 
transformation is not a symmetry of the theory, it does preserve some of its 
properties, and therefore can be used to identify phenomena on the black 
hole side with phenomena in the gauge theory. 



C.l Boosted black string 

We start with the black string solution in eleven dimensions 



r, 



7 



J.7 



h{r) = l-^. (51) 
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If the eleventh dimension z is compactified on a circle of radius R one can 
reduce this metric to get the ten dimensional metric of a Schwarzschild black 
hole in type IIA string theory. The black hole mass and entropy are given 
by (/p is the eleven dimensional Planck length) 



M 



^sRrl 



S ^ ^ (52, 



where fls — ^'^ 



r(9/2) • 
Following [O, IT^I we boost in the z direction (in the covering space) with 



boost parameter a, and recompactify the new z direction on a circle with 
radius 

R' = -^. (53) 

cosh a 

This radius is chosen so that the resulting ten dimensional charged black hole 

will have the same entropy (^) as the Schwarzschild black hole. 

In string frame the charged black hole solution is 

ds^ = -r^'\r)h{r)dt' + f'\r)[h-\r)dr^ + r^d^l] 

, , ^ rl sinh a cosh a , , 

r ' + Tq smh ot 

We will eventually take the limit a — >■ oo with tq fixed (this limit was also 
considered in |51|, |5^, |5^), so we won't distinguish between sinh a and cosh a. 
In this limit the number of unit charge DO-branes is given by dividing the 
ADM mass of the black hole by 1/i?', 

A^ = -MR. (55) 



Now let's consider the string theory parameters that come from this com- 
pactification. Our conventions are 

R = 9sls ^' = | (2%7)l/3 - (56) 
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One finds (string parameters are labeled by a tilde) 
~2 Ip cosh a 



i?(2%7)l/3 

Ip cosh a / 

9ym = ^gJs = ^^2/6 ^Q^h «• 

The string length grows with a, so we rescale all lengths by 1/ cosh a. 
This gives new parameters (labeled by ') 

fA2 ^ V cosh-V 

g's = 9s 

[9ym) - -^ 4^2 ^^^) 



J 



r/ cosh a. 



Now taking the limit a ^ oo brings us into the gauge theory regime with a 
finite value for the Higgs vev 



J 



U^77^,= r/ll (59) 



(I's 

On the supergravity side the change of scale means that 



ds'^ = ^ds^ . (60) 

cosh a 



Using equations (^,0,^,^) we can write 
rl sinh^ a _ d^Ng'^^ 



r7 f/7/;4 

.7 n tt7 



^ = ^«of/-^M,;t. = ^ ^ KU) - 1 (61) 



where do = T{7/2)2'^7i^^^ and ao 



r(9/2)2"7r"/^ 



As a — > oo note that l^ — > 0, so the metric becomes (defining t = 
1 1 cosh a, e^ = Qym^) 

TT7/2 jl/2 

ds' = lJ[-t__hiU)dt' + h-\U)%^dU' + d'J'eU-'/'dnl] . (62) 
do e ^ 

This is indeed the near horizon geometry of a 0-brane black hole. This 



analysis is closely related to [Q 
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C.2 Parameter identification 

To give a precise identification between the Schwarzschild black hole and 
the gauge theory we recall that a near extremal black hole has the following 
energy density, temperature and entropy density [El]- 



(63) 



Here e^ = Qym^ ) Uq is the position of the horizon in Poincare coordinates, 
and p + 1 is the worldvolume dimension of the gauge theory. The constants 
appearing in these equations are 



E = 


- a;'e-'N^Ur 


eT = 


- CpUQ 


S = 


= h,e''N'ut'^/' 



dp 



'9-Pao11-2« H^ 



9-p 

c, = (7-p)^/^r-i(^)2-V^ 



K = (7-j))-^/^r-^(^)2^-V^ (64) 



The above equations are only reliable if the supergravity approximation 
is valid near the horizon of the black hole - that is, if the dilaton is small 
and the curvature is small in string units. The condition for small curvature 
near the horizon is Uq'^ <^ e^, and the point where the curvature is order 
one is T ~ [/q. As a passing observation note that while the thermodynamics 
depends on p, one has a relation E ~ e~^A^^T^f/g which holds independent 
of p, aside from the exact numerical factor. 

Now we can make identifications between the Schwarzschild black hole 
and the gauge theory. We consider the case p = 0. On the Schwarzschild 
black hole side there are three parameters M, R, Ip. On the gauge theory 
side there is N, E, Qym- By equating the thermodynamics we obtain the 
relations 

9 

E = —M 
16 

2 _^2Srf^R^ 
^^*^ " 4vr2 II 

N = -RM 
8 

7 

-'gauge -'Schwarzschild ~. \ '-''-*/ 

Anrn 
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Uo — — 






>^p"a,iip^p >^I 



gauge ^ Schwarzschild 



M = 


-- E = N^e-^/'T^^/' 


R = 


= A^-ig^/ST-i^/s 


S = 


= A^2g-6/5j.9/5 


11 - 


= N-V^e^/^5rjn-U/5 



C.3 Transition points 

Let's see what the different regimes of the gauge theory correspond to for the 
Schwarzschild black hole. It is convenient to use (we drop numerical factors) 



(66) 



There is one issue we would like to comment on. The specific heat of 
a Schwarzschild black hole is negative, while the gauge theory clearly has a 
positive specific heat. Nonetheless their thermodynamics can be identified, 
since we hold different quantities fixed when we take derivatives with respect 
to the temperature (the temperature is the same on both sides) . In the gauge 
theory we hold A^ and Qym fixed, but on the supergravity side we hold R and 
Ip fixed. This is responsible for the difference. Thus the negative specific heat 
of the Schwarzschild black hole is not in contradiction with its description in 
terms of thermodynamics of an ordinary field theory (see also ||1CI|| ). 

Now let's look at some transition points: 

1. When R = tq the black string becomes unstable and turns into an 
eleven dimensional black hole |^. Our formulas are only appropriate 



for R < tq. At this point one finds 

Si = N 



R, = iv5/9e-2/3 

Ti = e2/3Ar-5/9 (67) 

/2 = e-^/3^2/9_ 



p 



This is the BFKS point ||, |T5|. 



2. At i? = Ip the black string is better described as a ten dimensional 
black hole and the string theory has Qs = 1- At this point 

^2 = Ar8/7 

i?2 = Ari/3e-2/3 (68) 

T2 = e2/=^Ar-i°/2i _ 
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3. Consider the point where the Schwarzschild black hole has a horizon 
size of the order the string length, rg = /^ = l^/ R. We find the corre- 
sponding gauge theory parameters 

^3 = N^ 

T3 = e^/2 (69) 

1 

9s - j^- 

This is the point where the curvature at the horizon of the charged 
black hole becomes of order one. Beyond this point the effective cou- 
pling of the gauge theory is small. We have argued in this paper that 
at this point there is a phase transition. We now see that this tran- 
sition occurs at the Horowitz-Polchinski correspondence point, where 



the Schwarzschild black hole turns into an elementary string state |47| 
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